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Abstract — A novel class of auction-based games is formulated 
to study coordination problems arising from charging a popula- 
tion of electric vehicles (EVs) over a finite horizon. To compete for 
energy allocation over the horizon, each individual EV submits 
a multidimensional bid, with the dimension equal to two times 
the number of time-steps in the horizon. Use of the progressive 
second price (PSP) auction mechanism ensures that incentive 
compatibility holds for the auction games. However, due to the 
cross elasticity of EYs over the charging horizon, the marginal 
valuation of an individual EV at a particular time is determined 
by both the demand at that time and the total demand over the 
entire horizon. This difficulty is addressed by partitioning the 
allowable set of bid profiles based on the total desired energy 
over the entire horizon. It is shown that the efficient bid profile 
over the charging horizon is a Nash equilibrium of the underlying 
auction game. An update mechanism for the auction game is 
designed. A numerical example demonstrates that the auction 
process converges to an efficient Nash equilibrium. The auction- 
based charging coordination scheme is adapted to a receding 
horizon formulation to account for disturbances and forecast 
uncertainty. 

Index Terms — Electric vehicles; progressive second price auc- 
tion; incentive compatibility; cross elasticity; game theory; Nash 
equilibrium. 


I. Introduction 

Vehicles that connect to the electricity grid to recharge, 
referred to generically as electric vehicles (EVs), offer a 
range of potential benefits, including reductions in reliance on 
liquid fuels and in pollutant emissions, and increased energy 
efficiency [ 1 ]— [3] . It is therefore anticipated that EV sales will 
substantially increase over the next few years [4]. If such 
growth does eventuate, it will become necessary to account 
for EV charging patterns in grid operation, as argued in [4]- 
[7] and references therein. 

Accommodating large numbers of vehicles on the grid will 
require coordination of EV charging so that their power and 
energy requirements can be optimally and robustly satisfied. 
This is a challenging control problem. Work on analyzing 
EV charging schedules, and their effect on utilities, began 
in the 1980s [8] . Recent work is extensive and includes [9]— 
[11] which formulate EV charging control as constrained op- 
timization problems, and [12] which develops an EV dispatch 
algorithm in the context of a day-ahead electricity market. 
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Centralized coordination faces numerous difficulties, from 
computational complexity to the loss of EV decision-making 
autonomy. Many distributed coordination methods have been 
proposed to address those challenges, including [11], [13]- 
[22] and references therein. In a general sense, that work is 
structured around individual players determining their optimal 
charging strategy over the charging horizon with respect to 
either the total demand of the other players, or the system 
(clearing) price, which is based on the total system demand. 
More specifically, a hierarchical structure is considered in 
[14], [15], [20], [21] for scheduling EV charging. Each EV 
determines its preliminary charging (load) profile by solving 
an individual optimization problem with respect to the latest 
forecast of the system clearing price. The clearing price is 
then updated to take into account the latest charging profiles 
of the individual EVs. For the update process designed in 
[14], the resulting strategies asymptotically approach a Nash 
equilibrium (NE) as the EV population increases to infinity. 
The resulting NE is nearly socially optimal. The market design 
in [21] considers both energy and reserve capacity 1 , and also 
incorporates distribution network losses. It is shown that in 
such a market setting, the NE may not coincide with the 
socially optimal solution. To establish a tractable formulation, 
only energy markets will be considered in the remainder of 
the paper. Also, network losses are not modelled. 

Most of the distributed methods cited above are quite 
distinct from the economic generation dispatch that under- 
pins deregulated day-ahead electricity markets [23]. To eco- 
nomically dispatch generation, auction mechanisms, such as 
uniform market-clearing-price [24] and pay-as-bid [25], have 
been widely adopted in electricity markets around the world 
[26]. Each generating unit submits to the ISO their bids 
over the forward market period (typically 24 hours), with 
bids consisting of pairings of minimal selling price and 
maximum supplied electricity for each market subinterval. 
The ISO dispatches the generation requirements among units 
based on their submitted bid profiles. However these auction 
mechanisms do not achieve incentive compatibility and usually 
cannot attain the efficient solution [27]. In contrast, this 
paper studies EV charging coordination over multiple time 
intervals under an incentive compatibility mechanism [28], 
[29]. The paper utilizes a progressive second price (PSP) 
auction mechanism, designed by Lazar and Semret [30], [31] 
and initially applied in the allocation of network resources. 

In a single divisible resource allocation problem under the 
PSP auction mechanism, each player only reports a two- 
dimensional bid. This bid is composed of a maximum amount 
of demand and an associated buying price, and is used to 

'Reserves refer to generator and/or load capacity that is available to 
compensate for sudden changes in energy supply or demand. 
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replace the player’s complete (private) utility function. Under 
the PSP mechanism, the money transfer (or payment) of a 
player measures the externality that they impose on the system 
through their participation. This concept will be formalized 
in Section III-B. As analyzed in [30], [31], the PSP auction 
mechanism is a VCG-style auction [32]— [34]. Therefore incen- 
tive compatibility holds, ensuring that all players submit truth- 
telling bids, and resources are allocated efficiently. Under this 
mechanism, as verified in [31], [35] in the context of single- 
unit network resource allocations, the efficient bid profile is a 
NE. 

In formulating their bids, players must consider tradeoffs 
between energy costs that vary over the charging horizon, the 
benefit derived from the total acquired energy, and battery 
degradation. Individual EVs are therefore inter-temporal cross- 
elastic 2 loads, as defined in [36]. This results in an auction- 
based allocation of a collection of divisible resources, where 
electric energy at each time-step of the horizon is a separate 
divisible resource. Consequently, each EV must submit a bid 
that has dimension double the number of divisible resources 
to be shared (equivalently double the number of time-steps in 
the charging horizon). Such auctions have received limited 
attention in the literature. It will be shown that a player’s 
marginal valuation for electric energy at a particular time is 
dependent upon both the amount of energy requested at that 
time and the total energy request over the entire charging 
horizon. A key contribution of the paper is to show that 
the efficient set of EV bids over the charging horizon is a 
NE of the underlying auction game. However, due to the 
cross elasticity of a bid over the multi-step time horizon, it 
is infeasible to directly verify this NE property using analysis 
that is applicable for a single-resource auction game [31], [35]. 
An alternative approach is proposed in this paper. 

In order to address cross elasticity, the set of bids of a 
player is partitioned into a collection of subsets, each of 
which is composed of bids that possess the same total desired 
electric energy over the horizon. Consequently, cross elasticity 
is eliminated for bids within each subset. For such bids, the 
marginal valuation at each time-step includes a variable part 
determined by the amount of energy requested at that time and 
a. fixed part that is identical for all bids in that subset. With this 
construction, the NE property of the efficient solution can be 
established by verifying for each subset whether any player 
can benefit by unilaterally deviating from their efficient bid 
profile. 

The paper is organized as follows. Section II establishes 
the problem structure by formulating a class of EV charg- 
ing coordination problems over a multiple time-step horizon. 
Distributed charging under the PSP auction mechanism is 
introduced in Section III. Section IV shows that the efficient 
(truth-telling) bid profile is a NE of the underlying PSP 
auction game. In Section V, an update process is designed to 
implement the PSP auction, and an example illustrates that this 
process converges to the efficient NE. Section VI concludes 
the paper and provides a discussion of future work. 

2 Cross elasticity refers to the ability of EVs to move charging demand from 

one time interval to another. 


II. Electric Vehicle Charging Coordination 
Formulation 


A. Charging and cost models 

This study focuses on coordinating the charging of a 
population of EVs, denoted by A f, finite charging horizon 
T = {foi • • ■ , to + T — 1}, with to representing the initial 
time. For each EV, n £ Af, the energy delivered 3 over the 
t-th time interval is denoted by x nt , and the battery state of 
charge (SoC) evolves according to: 

Sn,t + 1 = Snt 4 7T Xnt-, (1) 

t'n 

where (-)„ is the battery capacity and s nt is the normalized SoC 
for the n-th EV at time t. An admissible charging strategy , 
x n = ( x n t,t £ T), satisfies the constraints: 



when t gT u 
otherwise 


with ^2 x nt < r„, (2) 

ter 


where T n C T denotes the charging horizon of the ri-th 
EV, r n = 0„(s n — Sno) gives the maximum energy that it 
can receive, and 0 < s n o < s n < 1 give the normalized 
minimum (initial) and maximum SoC, respectively. The values 
for T n , r n and s n o follow from the driving style and vehicle 
battery capacity, see for example [37]. The set of all possible 
admissible charging strategies is denoted by X n . Also, define 
the collection of admissible charging strategies for all EVs by 
x = {x n \ n £ Af), with its corresponding set being X. 

As specified in (2), the charging demand of each EV, n £ Af, 
at each time interval is elastic, but charging over the time 
horizon is coupled through the maximum energy F n that 
can be delivered during that time. As a consequence, EV 
coordination is inherently a problem of scheduling demand 
that exhibits inter-temporal cross-elasticity [36]. 

The utility function of the n-th EV, for a charging strategy 
x n , is given by: 

tXn (®n) — 'y ' fn {Xnt) ^ ^ ' X n t , (3) 

teT teT 


where S n > 0 is a fixed parameter, and /„(•) denotes the 
battery degradation cost of the n-th EV. 


Remarks: 

• This work considers LiFeP04 lithium-ion batteries which 
have been widely used in a variety of EVs. Key char- 
acteristics of this type of battery, including state of 
health, growth of resistance, and cycle life, are effected 
by charging behavior over many cycles. A degradation 
cost model for this type of battery cell is formulated in 
[38], based on the evolution of battery cell characteristics 
developed in [39], [40], This degradation model expresses 
the energy capacity loss per second of a cell with respect 
to the charging current and voltage. Using this model, the 
relationship between degradation cost and charging rate 
x n t can be developed. 


~’It is assumed that the energy is delivered at a constant rate (power) over 
each time interval, and that the time intervals are of unit length. Therefore 
the charging rate is also given by x n t- 
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• The second term in (3) captures the penalty cost due to 
not fully charging the EV over the time horizon, with 
S n weighting the relative importance of delivering the 
maximum energy during charging [41]. 

• The utility function w n { x n ) therefore establishes the 
tradeoff between the battery degradation cost and the 
benefit derived from delivering the full charge. 

For a collection of admissible charging strategies x , the 
system cost is given by: 

Js(x) = 'Y^CtiyDt + ^2 Xnt) ~ ^2 W n( x n), (4) 
teT nGM nGM 

where c*(-), D t and D t + zLnc.v x nt denote the generation 
cost, the aggregate inelastic background demand and the total 
demand at time t respectively. It is assumed that a forecast 
for D t is available over the charging horizon. For the example 
presented in Section V-C, the demand profile D = (D t ; t G T) 
is given by the aggregate demand for a typical summer day in 
the Midwest ISO region of North America. 

The system cost function (4) considers tradeoffs between 
total generation cost, aggregate battery costs and the penalty 
for deviating from full charging. This contrasts with the gen- 
eral literature, for example [14], [15], [18], where the objective 
is to achieve valley-filling. While valley filling minimizes the 
total generation cost, it may result in high battery degradation 
costs across the EV population. 

B. Efficient charging 

It is desirable to determine the collection of efficient (so- 
cially optimal) charging strategies x ** that minimizes the 
system cost (4). This centralized EV charging coordination 
problem can be formulated as the following optimization 
problem: 

Problem 1: 

min J s (x ), 

such that x satisfies constraints (2) for all n £ TV. ■ 

The efficient charging strategy x ** of Problem 1 can be 
characterized by its associated KKT conditions. Firstly, the 
Lagrangian can be written: 


(A 1 ) ct(y) is monotonically increasing, strictly convex and 
differentiable on y\ 

(A2) fn{x), for all n £ TV", is monotonically increasing, 
strictly convex and differentiable on x. 

Remarks: 

• The generation cost ct(-) is widely assumed to be a 
convex function of total generation, see for example [42]— 
[44], 

• The battery degradation cost /„(•) is governed by the 
chemical processes inherent in charging. It is shown in 
Fig. 7 of [39] that growth of battery resistance, hence the 
fade of battery energy capacity, is generally increasing 
and convex with respect to charging rate. This provides 
some justification for (A2) since f n measures the cost 
related to the fade of battery capacity with respect to 
charging rate. 

Lemma 2.1: The collection of efficient charging strategies 
x** for Problem 1 is unique. 

Proof. Under Assumptions (Al, A2), the cost function 
J s (x) is strictly convex and differentiable. Also, the con- 
straints (2) determine a convex domain. Therefore Problem 1 
is a strictly convex optimization problem. Thus there exists a 
unique solution. ■ 

This centralized charging coordination strategy can only 
be effectively implemented when the system has complete 
information and can directly schedule the behavior of all 
EVs. In practice, however, individuals are often unwilling 
to share their private information with others. Furthermore, 
transmission of complete information may incur excessive 
communications, and centralized control might be computa- 
tionally infeasible. Thus this paper focuses on the development 
of a distributed control method that is based on the progres- 
sive second price (PSP) auction mechanism, which has been 
applied in [30], [31], [35] for efficient allocation of single- 
unit network resources. However, because EVs are scheduled 
over a multiple time-step horizon, charging coordination is a 
multi-unit resource allocation problem. 

III. Distributed EV Charging Coordination under 
a PSP Auction Mechanism 


L ( x , \ ) — J s ( x^j T ^ ' X n ^ ^ ' x n t , 

neV teT 


where X n is the Lagrangian multiplier associated with the 
constraint 'f2teT Xnt — V' f rom (2). The KKT conditions 
for Problem 1 are therefore given by: 


-L(x, A) > 0, x nt > 0, 


-L(x, X)Xnt = 0, (5a) 


dXnt ’ ' ' dXnt 

^x nt -r„<0, A„ > 0, A„^^a:„t - Tnj = 0, (5b) 

ter ter 

for all t £ T and n £ TV, where: 

d q 

fff d/{x, A) — (,-Dt T 'y ' x n t) fff 4* An • (5c) 

nt neN nt 


A. Bid profiles of individual players 

Each EV, n £ TV, submits a 2T-dimensional bid, b n = 
( b nt ,t £ T), where: 


bnt ={/3nt,dnt), with 


dnt 



t £ Tn 

otherwise 


and ^2 dnt < r„ 


teT 


specifies the price f3 nt that player n is willing to pay for 
energy at time t and the maximum electrical energy d n t that 
is desired at that time. The corresponding feasible allocation 
x n = ( x n t,t £ T) with respect to b n must satisfy: 


0 < x nt < d nt , Vt S T. (6) 


Assumptions: The following conditions apply throughout the Let B n denote the allowable set of bids for player n, so that 
remainder of the paper: b n £ 
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Each player’s revealed utility function, denoted by 
w n (x n (b n ); b n ), is defined as: 

' l Xn{x n (b n ); b n ) = ^ ^ Pnt%nt- 

ter 

The revealed system cost with respect to a bid profile b = 
( b n ,n £ A f) is given by: 

J )x[b) . b ) — ^ ^ Ct (D t “t - ^ ^ %nt ) ^ (a? n )b n ) 1 . 

teT nGW" neM 

(7) 

Auction-based EV charging allocation can be written as the 
following optimization problem: 

Problem 2: 


J*(b)= min J(x{b);b). 

Constraint (6) 

The objective of the auctioneer is to assign an optimal allo- 
cation x* (b) with respect to a bid profile b to minimize the 
revealed system cost given by J. ■ 

Unlike the single-sided auctions considered in [30], [31], the 
cross elasticity inherent in Problem 2 suggests that the optimal 
allocation x* (b) over the charging horizon T depends upon 
both the bid profile b over the entire horizon and the generation 
cost Cj(-). However, the following lemma shows that the 
optimal allocation x* at time t, is completely determined by 
the bid profile b t and Ct ( ■ ) at only that time. 

Lemma 3.1: Suppose x*(b) = (a :^,t £ T) is the optimal 
allocation subject to bid profile b. Then: 

x 1(b) = Xt(b t ), for all t £ T. 


Proof. For notational simplicity, consider: 


S(b t ) = {x t ; s.t. 0 < x nt < d nt , for all n £ AT}, 

bt {Xf. bfj = Ct(D t -f ''y ' Xnt) ^ ^ PntXnt- 
neJV n£ M 

It follows from (7) that: 

min J(x(b);b) = min ht(x t :b t ) 

Constraint (6) V x t eS(b t ) ,teT V 

= V min h t (x t ;b t ), 
tTr XteS{bt) 

where the last equality holds because the summation is sepa- 
rable in terms of x t . The desired conclusion follows. ■ 
The optimal charging allocation x*(b) of Problem 2 can be 
characterized by the associated KKT conditions. Firstly, the 
Lagrangian can be written: 

L a (x, cr; b) = J{x(b);b) + ^ - d nt ), 

neM tST 


where a n t is the Lagrangian multiplier associated with each 
constraint in (6). The KKT conditions for Problem 2 are given 
by: 


d d 

L a (x, cr; b) > 0, x nt > 0, -^—L a (x, cr; b)x nt = 0, 

(8a) 

( x n t d n t)cr nt — 0, 

(8b) 


\ ? 1 / w 7 /it v 

UXnt UXnt 

%nt d n t — 0? ®nt ^ 0? 


for all t G T and n G A f, where: 
d 

L a (x , cr; b ) = c' t (D t + ^ x kt ) - /3 nt + o nt . (8c) 
keM 


dx r 


It is now possible to establish a connection between the 
optimal charging strategies given by Problems 1 and 2. 
Lemma 3.2: Consider a collection of bids, 


Kt = = 


d 


8x r 


■Wn[X n ),X nt 


(9) 


for all n £ Af and t £ T. Then, under Assumptions (Al, A2), 
x*(b*) = x**. Also, 


Pit 


= c't(Dt + J2keM d lt ). if <t > 0 
< cj(A + E fce .V d fct)> = 


( 10 ) 


for all n £ Af and t £ T. 

This lemma is essentially the so-called fundamental the- 
orem of welfare economics [45]. Verification is provided in 
Appendix A. ■ 


Remark: Lemma 3.2 establishes two important properties: 

(i) It specifies a bid profile b * , given by (9), under which the 
optimal charging allocation x* of Problem 2 is efficient. 

(ii) At bid profile b* , EVs with an allocation larger than 
zero share the same marginal price as generation, which 
is larger than or equal to the price of EVs with zero 
allocation. 

Incentive compatibility holds under the PSP auction mech- 
anism [30], [31]. Therefore, a bid with price satisfying (3 n t = 
-g^—w n {d n ), for all t £ T, as is the case in (9), is the best 
choice among all possible bids. It follows from (3) that the 
truth- telling bid of the n-th EV is given by: 

Pntidnt',^2 d ns ) = -fn(d nt ) +2S n (T n ~ ^ d ns ^ . (11) 

sGT sGT 

This implies that an EV’s marginal valuation at each time-step 
is determined by both its electrical energy request d nt at that 
time and its total energy request E seT dns over the entire 
multi-period charging horizon. 


B. Calculation of EV payment and payoff 

The payment incurred by each EV will be specified with 
respect to the allocation law defined by Problem 2. Each 
EV’s payment is exactly the externality imposed on the system 
through its participation in the auction. For the n-th EV, this 
is given by the system-wide utility when the n-th EV does 
not join the auction process, minus the system-wide utility 
(but excluding the contribution of the n-th EV itself) when 
the n-th EV joins the auction. 

To express this payment, it is convenient to introduce a 
slight abuse of notation by writing the collection of bids as 
b = ( b n , b- n ), where b- n = {b kl k £ Af\{n}). The payment 
of the n-th EV, for a bid profile b, is then given by: 

T n (b) = -J*(On,b- n )-^-r(b)-Y / p nt X* nt (b)j , (12) 
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where (0„, b_ n ) denotes the bid profile without the 77 -th EV’s 
participation, i.e., with the bid d n t replaced by d n t = 0 for all 
t G T, and x*(b) is the optimal charging allocation given by 
Problem 2, with respect to b. Thus by (7), ( 1 2) and Lemma 3.1, 
the payment of player n at time t can be expressed as: 

Tnt{b t ) = - C t (D t + ^ Zm7") + Ct(A + Y X kt( b t)) 
m^n k£j\f 

+ Y P mt ( X ™t n ~ X *mt( b t )) , 

m^n 

where x*fi n denotes the optimal charging allocation of EVs 
m ^ n given by Problem 2 with respect to (0„, 6_ n ) t . (Recall 
from Lemma 3.1 that allocations at time t are unrelated to 
other times.) The total payment of player n is given by: 

T n ( b ) = Y T nt( b t)- (13) 

ter 

The payoff function of the n-th EV is given by the difference 
between the EV’s utility and its payment: 

u n (b) = w n (x* n (b )) - T n (b). (14) 

This payoff function provides the basis for defining a Nash 
equilibrium for the PSP auction game. 

Definition 1: A collection of bid profiles if is a Nash 
equilibrium for Problem 2 if: 

u n (b^,b°_ n ) > u n (b n ,b°_ n ), 

for all b n G B n and for all n G A f. That is, no EV can benefit 
by unilaterally deviating from its bid profile 6° . ■ 

C. Related work on EV charging games 

1 ) Hierarchical EV charging games: Many of the dis- 
tributed schemes that have been proposed for scheduling EV 
charging have adopted a hierarchical structure, see for example 
[14], [15], [20], [21]. Each EV first determines its optimal 
charging schedule with respect to the expected energy price 
that is broadcast by the ISO. The price is then updated based 
on the latest charging schedules attained by the population of 
EVs. This process repeats until price updates become negli- 
gible. Typically, the outcome of such a scheme approaches a 
NE asymptotically as the size of the EV population approaches 
infinity [14]. In contrast, the auction game formulated in this 
paper achieves a NE for small collections of EVs. 

2 ) Efficiency of NE for auction games: The tradeoff be- 
tween efficiency and risk is considered in [22] in the context 
of scheduling load over a finite time horizon in oligopoly 
electricity markets. In non-cooperative schemes, the system 
may suffer a loss, known as the price of anarchy, which relates 
to the difference between the efficient solution and the NE. The 
PSP auction proposed in this paper circumvents that loss by 
ensuring the NE is efficient. 

Efficient games have been studied in a variety of related re- 
source allocation problems. The energy consumption schedul- 
ing game formulated in [18] does not employ user utility 
and considers a linear payment rule that differs from VCG. 
It achieves a globally optimal NE that minimizes the total 
generation cost. The efficiency of single-unit network resource 


allocation is studied in [30], [31], [35], though EV charging 
coordination is essentially a multi-unit resource allocation 
problem. 

IV. Efficiency of the Charging Coordination PSP 
Auction Game 

This section establishes that the efficient solution given by 
Problem 1 is a NE of the PSP auction game. Budget balance 
under the PSP auction is considered in Section IV-D. Finally, 
it is shown in Section IV-E that for a single-interval auction 
game, where T = 1, the efficient NE is unique and so there 
in no efficiency loss. 

Suppose b* is the bid profile specified in Lemma 3.2, such 
that the corresponding optimal charging allocation is efficient. 
It will be shown in this section that b* is a NE for the 
underlying auction game. By Definition 1, this implies: 

U n (b* n ,b*_ n )>U n (b n ,b*_ n ), (15) 

for all b n G B n and for all n G TV". 

Due to the cross elasticity arising from the second term 
in (11), directly verifying that the efficient bid profile is the 
best for every player is infeasible. Accordingly, the approach 
developed in [35] for auction games of a single divisible 
resource is not applicable. In order to overcome this difficulty, 
an alternative approach will be developed. This involves parti- 
tioning the set of bid profiles B n into a collection of subsets: 

Bn{A) = jf>„ G Bn, s.t. Y dnt = a|, (16) 

ter 

each of which is composed of admissible bid profiles that 
possess a common total A for the desired energy over the 
charging horizon 77 The set of bids is then given by: 

B n = U B n {A), 

Ae[ o,r„] 

noting that B n (A)f]B n (A) = 0 whenever A A. For all 
bids in a particular subset B n (A), it follows from (11) that the 
marginal valuation price /3 n t at any time t includes a variable 
part —f' n (dnt) that is dependent upon the request d n t at that 
time, and a fixed part 2 6 n (IV, — A) that is identical for all bid 
profiles in B n (A). Therefore, the cross elasticity over the time 
horizon is avoided. 

By Definition 1 and the specification of B n (A), it is suffi- 
cient to show that b* is a NE, if for every fixed A G [0, L,,]: 

u n {b* n , b*_ n ) > u n (b n ,b*_ J, for all b n G B n {A), (17) 

and for all n G TV". This is easier to verify than (15), since 
the difficulty associated with cross elasticity is avoided when 
considering each specific subset B n (A) with fixed A G [0, T n ]. 

It will be shown in Section IV-A that (17) holds when A > 
ff tr-j- ( K lt , and Section IV-B considers the case where A < 
fftr-j- df t - These results build on the following lemma. 

Lemma 4.1: Suppose /3 n t(d n t', A) is the bidding price given 
by (11) for player n at time t, but with the summation ff , d n t 
in the second term replaced by A. Then f3 nt (d n t ; A) satisfies 
the properties: 

Pnt(dn t ; A) > nt (d 2 nt ; A) > 0, with d} nt < d 2 nt , for all A, 

(18a) 
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Pnt{dnuAi ) > P n t(d nt ; A 2 ), with A x < A 2 , for all d nt . 

(18b) 

Proof. By (11), fi n t(d n f,A) = 26 n (T n -A)- f' n (d nt ). It 
is straightforward to verify (18) under Assumption (A2). ■ 

Lemma 4.1 implies that /3 n t(d n t , A) decreases with increas- 
ing d n t and A, as illustrated in Fig. 1. 



Fig. 1. Illustration of /3nt(d n t', A) with respect to d n t and A. 


A. Verification of (17) when A > YhteT d* t 

The first step in verifying (17) is to show that all EVs m £ 
M \ {n} are fully allocated when A > YteT d nt- This is 
established by Lemma 4.2. The main result then follows as 
Theorem 4.1. 

For notational simplicity, let x* and x denote the optimal 
allocations with respect to b* and (b rl . b*_ n ) respectively. 

Lemma 4.2: If A > ]P tg7 -d* t , t ^ len 

Xmt = d* m t, for all m £ fif \ {n}. 

Proof. If d^t = 0 for any m £ Af\{n}, then x m t = x* mt = 
0 and the desired result is obtained trivially for that m at time 
t. It will be assumed that > 0 in subsequent analysis. 

Case I, when d n t > d* t . From Lemma 4.1, since d n t > d] \ t , 
(18a) gives: 

Pnt(d n f, J2 <t) < Pnt(d* nt ',J2 d nt) = P*nt, 

teT teT 

and because A > Y2teT d nt, (18b) gives: 

(dnt = (dnt (dnt 5 A ) A (3 n t (d n t i ^ ^ ^nt)’ 

ter 

This implies (d nt < (d* t , as illustrated in Fig. 2, with equality 
holding only when d n t = d* t and A = YteT d nt- l n that 
special case, the bids (b n ,bf n ) t and b * coincide, so x m t = 
x^t = d mt as desired. The following analysis considers fd nt < 
!3* t . Given the earlier assumption that d* nt > 0 , Lemma 3.2 
indicates that (d* t < Hence, fi nt < (d* t < (d* mt . 



Fig. 2. Illustration of (3 n t and /3* t when A > 


If x nt > 0 then by (8a), g§-^L a (x, a; (b n , b*_ n )) = 0. 

Because a nt > 0, (8c) implies c' t (D t + J2k&jV^kt) < (d n t ■ It 
follows that, 

c' t (D t + Y, *fct) < Pnt < P*mt = c' t (D t + Y *«)> (19) 

fceAf fceAf 

and hence that JfkeM Xkt < YlkcAf x lt- Therefore the total 
energy allocated to all the EVs at time t decreases under 
the bid profile (6 n ,6L„). Furthermore, (19) gives c' t (D t + 
'YkcAT x kt) - fdmt < 0, so (8c) implies that o mt > 0. The 
complementarity condition (8b) then ensures x mt = d^t for 
all m £ M \ {n}, which is the desired result. 

If x n t = 0, then Xkt < d% t = x% t for all k £ J\f. Assume 
Xmt < d* mt for some m £ AC \ {n}. Then J2keAT x kt < 
Y2keN x tv and so 

c' t (D t + Y *w) < <4(A + Y Xk t) = P*mt- (20) 

keAf keAf 

Again, (8c) implies o m t > 0, which is inconsistent with the 
assumption x m t < d* rlt because of complementarity (8b). 
Therefore, x mt = d* , for all m £ J\f \ {n}, as desired. 

Case (II), when d n t < d* t . Because x n t < d n t < 

d*nt = Xnt, then 'EkGAfXkt < YktAf X kV and S0 ( 20 ) holds - 
It again follows from (8c) that o m t > 0, and hence the 
complementarity condition (8b) ensures that x m t = d * rlf , as 
desired. ■ 

The main result for this section can now be established. 
Theorem 4.1: Under Assumptions (Al, A2), (17) holds 
when A > J2teT d nt- 

Proof. Considering the n - th EV, the first step is to com- 
pute the difference between the payoff given by the optimal 
strategy u n (b*) and that obtained from an alternative strategy 
u n (b n ,b*_ n ). Using (14), these payoffs are given by, 

u n (b*) = w n (d* n ) - r n (b *), (21a) 

u n (b n , b *_ n ) = w n (x n ) - T n (b n ,b*_ n ). (21b) 

Hence, the difference is: 

A u n = u n (b*) - U n (b n ,b*_ n ) 

= w n (d* n ) - r n (b*) - (w n (x n ) - Tn(b n , bl n )J 

= W n (d* ) - W n (x n ) + (^J*(b n ,b*_ n ) + Y IdntXnt 

teT 

-J*(b*)-YPnt d nt)’ 

teT 

where the final equality follows from (12). Straightforward 
analysis, using the result from Lemma 4.2 that x m t = d* rlt for 
all m £ M \ {n}, gives: 

A U n — tlJ ri (d n ) W n (Xn) "(“ ^ -\~ 'y ( d mt A X n t) 

teT m^n 

-Ct(D t +Y d *kt)}- 

keAf 

Also, using Lemma 4.2, 

J s (x*) - J s (x) = J s (d*) - J s (x n ,d*_ n ) 
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— E {^{Dt + X] d kt) ~ c t{D t + E ^mt + x iit)} 

tGT k€.Af m^n 

— ^ ^ w k(d'fc) + ^ ^ Wm^dyn) 4 - W n (x n ) 
k£Af m^n 

= E { c *(-^ + E — c *(-^t + E ^* mt + x ™*)} 

teT /ceA/ - m^n 

^'n (tf n ) T 'Utn(x n ) 

= -Au n . 

Since a:* is the efficient allocation, J s (x*) < J s (x), so 
Au n = J s (at) - J s (x*) > 0. 

This implies that player n cannot benefit by unilaterally 
changing its bid 6* to any other bid b n £ B n (A) with 

A A EtGT d nt- ® 


B. Verification of (17) when 0 < A < EteT dl t 

This section considers a fixed A satisfying 0 < A < 
EteT d nt- Firstly, Lemma 4.3 establishes the n-th EV’s 
optimal bid in the set B n (A) and under the constraint. 


0 < d nt 



when d* nt > 0, 
otherwise, 


for all t £ r. (22) 


Theorem 4.2 then shows that the bid of Lemma 4.3 remains 
optimal when the constraint (22) is relaxed. Using this result. 
Theorem 4.3 finally establishes that (17) holds when 0 < A < 

Et£T d nf ^ ^ ^ ^ 

Lemma 4.3: Consider a bid 6* = b*(A) = ((/3* t , d* t ), t £ 
T), with A £ [0, EteT d nt)> suc h that 

K = argmax u n (b n ,b*_ n ). (23) 

£n£Z3„(A) 

Constraint (22) 

Let x* = (scl t , k £ A/", t £ T) denote the optimal allocations 
with respect to (t>l,b*_ n ). Then, 

= <C = <C for all 771 £ TV \ {n}. (24) 


Furthermore, define the function, 

gnt(d) A Ct{D t + E d mi + d) + fn{d). (25) 

m^n 

Then, under Assumptions (Al, A2), 6* satisfies the property: 


9nt(d* nt ) 



when d* nt > 0, 
when d* t = 0, 


for all t £ T, (26) 


where fi is a constant. ■ 

The proof of Lemma 4.3 is given in Appendix B. It will now 
be shown that the bid b* established in this lemma remains 
optimal when the constraint (22) is relaxed. 

Theorem 4.2: Suppose that b* is the optimal bid from 
Lemma 4.3. Then, under Assumptions (Al, A2) and with 
A € [0,E teT d nt)’ K satisfies: 

b* n = argmax u n (b n ,b*_ n ). 
bn&B n (A) 


Proof. Assume that A £ [0, EteT d lt)- Consider a bid 
bn = (( Pnt,d nt ),t £ T) £ B n (A) such that d nt > d* nt 

for some t £ T. Then the desired result follows if it can be 
proven that b n cannot be the optimal bid profile in the subset 
B n (A). 

Firstly, the following points hold: 

(i) P nt(d„ t ;A) > /3nt(d„t;'EteT d nt) because A < 
EteT dlt’ by (18b). 

(ii) c' t (D t + Em^n d *mt + d n t) increases with d n t under 
Assumption (AT). 

(iii) (/3* t ,d* t ) is the point at which Put {d nt ■, EteT d lt) and 
c t(A + E m^n d *nt + d nt) coincide, by Lemma 3.2. 

Using (i)-(iii) together with (18a) gives, 

0 <d* nt < fa, p* nt < Pl t , 

where {Pl t ,dl t ) denotes the point of intersection between 
fint {dnt 7 A) and c' t ( D t + Emjtn d ht +dnt)- This relationship 
is illustrated in Fig. 3. 



Fig. 3. An illustration of partitioned subsets 7 Zi, i = 0, • • • ,3. 

_ Let (Pliant) denote the bid with respect to A such that 
Pit = Pit > in other words Pl t = P* nt = p nt {Sl t , A). It follows 
from (18a) that dl t < Based on the ordering 0 < d* t < 
d\ t < d‘f lt , the set [0,oo) can be partitioned into four disjoint 
regions, 

^o = [0 ,dlt), 7 h = (dl t ,dlt\, 

n 2 = (dlt, dlt], fts = (&,+«>), 

together with the point d* nt . This partitioning is shown in 
Fig. 3, and is used in Appendix C. 

Suppose there exists a b n £ B n (A) such that at time t 2 , 
dnt 2 = d* nt2 . Then, because EteT d «* = E ter dlt = A < 
EXo 1 d lv and d^ < dl t f° r all t £ T due to (22), there 
must exist another time t \ such that d n ^ < d* t . Consider 
two bid profiles b n ,b n £ B n {A) such that, 


•+J 

V 

e 

'TS 

V 

. e 
■ 'TS 

(27a) 

d n t 2 ^ d n t 2 ^ d n t 2 — d nt2 , 

(27b) 

dnti H - d n t 2 — d n ti H - d n t 2 , 

(27c) 

dnt — dnti for all t ^ 

(27d) 

Then, it is shown in Appendix C that, 

^-n) ^ —n) ’ 

(28) 


Hence, such a b n cannot be the optimal bid with respect to 


Similarly, suppose there exists a b„ G B n {A ) such that at 
time t 2 , d n t 2 G IZ 1 UIZ 2 UIZ 3 . Following the same argument as 
previously, there must exist a time t\ such that d ntl < d* t . 
Again, consider two bid profiles b ni b n G B n ( A ) such that 
(27a),(27c),(27d) are satisfied, but with (27b) replaced by, 

dnt 2 ^ d n t 2 G d n t 2 ■ (29) 

Appendix C shows that (28) again holds, so b n cannot be the 
optimal bid with respect to b*_ n . This is the desired result. 

Moreover, because Assumptions (Al, A2) ensure differen- 
tiability of ct(-) and /„(.) for all n, the payoff function of 
every EV is continuous in the bid profiles. Therefore, it follows 
from (28) that the optimal bid for the n-th EV must lie in 
region 1Z Q for all f G T. ■ 

Theorem 4.3: Under Assumptions (Al, A2), (17) holds 
when 0 < A < EteT d nt- _ 

Proof. By Theorem 4.2, b* is the optimal bid with respect 
to bf n in B n (A), when A G [0, ffteT d nt)- Therefore, the 
best possible payoff for the n-th EV is, 

u n (K, b*_ n ) = w n (d* n ) - T n (b* n ,b*_ n ). (30) 

Using (24) and recalling that x* = d* gives, 

J s {x*) - J s {x*) = J s (d*) - J s (d* n ,d *_ J 

= E { c + E d *kt) - ctiPt +EC + E)} 

t£7~ /cGA/" m^n 

- E W k( d k) + E W m( d *m) + W»(<) 

k£Af m^n 

= E { C ‘(^ + E d kt) ~ CtiPt +J2 d *mt+ <t)} 

k£j\f m^n 

- ™n«) + Wn( d n) 

= - (u n (b*) ~ U n (6* , 6I„)^ = - A U n , 

where the final equality makes use of (21a) and (30). By 
Lemma 3.2, x* = x*(b*) is the efficient allocation solution, 
so J s (x*) < J s (x*). Therefore, 

A U n = J s (x*) - J s (x*) > 0, 

which implies that player n cannot benefit by unilaterally 
changing its bid b* n to any other bid 6* G B n (A) with 
A G [0,J2ter d nt)- ■ 

C. Existence of efficient Nash equilibrium 

Applying the analysis and results from Sections IV-A 
and IV-B, and under Assumptions (Al, A2), it can now be 
shown that the efficient bid profile b* specified in (9) is a Nash 
equilibrium of the EV charging coordination game under the 
PSP auction mechanism. 

Corollary 4.1: Under Assumptions (Al, A2), the efficient 
bid profile b* = (6* ; n G Af) specified in (9) satisfies the 
property: 

u n {b * n , b*_ n ) > u n {b n ,b*_ n ), for all b n G B n . (31) 

Proof. It was shown in Theorems 4.1 and 4.3, under 
Assumptions (Al, A2), that 

u n {b* n ,b*_ n ) > u n (b ni b*_ n ), for all b n G B n (A), 


holds when A > Ylt<£T d nt an ^ when 0 < A < J2teT d nv 
respectively. The desired result (31) holds since B n = 

U.4e[o,r„] &n(A). ■ 

D. Analysis of budget balance at the efficient NE 

As studied in [35], [46], [47], budget balance may not hold 
under VCG-style auction mechanisms. Define the surplus (net 
income) for a bid profile b as: 

m = E {^(b)}-E{ c t(A+E^t( b ))- c *( D ‘)}, 

neV tgr fee N 

which captures the aggregated EV payments minus the extra 
generation costs arising as a consequence of the EV popula- 
tion. Therefore, tp(b*) may not equal zero at the efficient NE 
b*. The sign of tp(b*) is established by the following lemma: 
Lemma 4.4: 

x!t{b*) > 0. 

The proof is given in Appendix D. ■ 

Lemma 4.4 establishes that the budget balance for b* will 
always be in surplus (positive). This surplus tp(b*) can be used 
to reimburse the services provided by the system operator. 

E. Efficiency loss of single -interval auction games 

Given that the efficient bid profile b* is a NE, efficiency 
loss can only occur if there exist multiple NE. The following 
theorem shows that for single-interval auction games, i.e., 
an auction game with T = 1, the efficient NE is unique. 
Hence, for such games, there is no possibility of efficiency 
loss. In this case, the efficient bid profile degenerates to 
b* = ((p*,d*),neAf). 

Theorem 4.4: b* = ((/?*, d*), n G A f) is the unique NE for 
PSP auction games where T = 1. 

Proof. Corollary 4.1 establishes that b* is a NE. Verifying 
Theorem 4.4 is equivalent to showing that every inefficient bid 
profile b cannot be a NE. This is undertaken in Appendix E. 

■ 

Ongoing research is considering efficiency loss for EV 
charging auction games over a multiple time-step horizon. 

V. PSP Auction Process for EV Charging 

Section IV established the existence of efficient Nash 
equilibria, under appropriate conditions, for the underlying 
EV charging coordination game. This section develops a 
bid-profile update process which motivates an algorithm for 
determining efficient NE. 

A process for determining an EV’s best bid, given the col- 
lection of bids for the other EVs, is presented in Section V- A. 
It will be shown that this process can be formulated as a 
dynamic programming problem. Section V-B then provides an 
algorithmic description of the update mechanism that governs 
the underlying auction game. A numerical example is provided 
in Section V-C to demonstrate this update mechanism and 
show that it achieves an efficient NE. 
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A. An EV’s best bid with respect to other EVs 

Recall that denotes the best bid of the n-th EV 

with respect to the bid profiles b_ n of all the other EVs, 

K(b- n ) = argmax u n (b n , b_ n ), 

where u n {b n ,b- n ) is the individual payoff of the n-th EV, 
as established in (14). However, due to the cross-temporal 
coupling arising from the summation term l2teT d nt in truth- 
telling bids, as identified in (11), it is impractical to directly 
determine the best response 6* (6_„) that is incentive compat- 
ible. This can be addressed by finding the best response when 
bids are constrained to possess a common total desired demand 
A = 'Zt^rdnt, and then optimizing over A £ [0,r n ]. The 
resulting optimization is given by, 

u n (b*,b- n ) = max max u n (b n ,b- n ). 
Ae[o,r„] b n &B n (A) 

The remainder of this section describes a dynamic program- 
ming approach to solving the inner optimal bidding problem 
that arises for each fixed total demand request A £ [0,r n ]. 

The dynamics associated with the physical charging process 
(1) can be rewritten, 

Sn,t-\- 1 {bn j it—n) — S n t {bn > b— n) T Xnt {bn , b__ n ) , (32) 

with t £ T, and where x n t{b n ,b-n) denotes the allocated 
charging rate of the n-th EV at time t, with respect to the bid 
profile 

Lemma 5.1: 

Consider a bid b n £ B n (A), for any fixed A £ [0, r n ]. Then 
the payoff function of the n-th EV has the summation form: 

tlnipni b~n ) = ^ ' Oj n {bt ) 5 n Q n (s n S n ) j (33) 
ter 

where a n {b t ) = -f n {Xnt(b t ))-T nt {b t ) and T nt (b t ) is defined 
in (13). 

Proof. From (14), 


u n (b) = w n {x* n {b n )) - r n {b) 

= Y { ~ fn(Xnt{bt)) - r nt (b t )| 

ter L 

^ 71 ( ^ ^ Xnt {bt ) , 

teT 


where equality holds by (3), Lemma 3.1, and (13). Then (33) 
follows directly from (32). ■ 

Let 7* = { f . .... I' — 1}, and define the value function 

V n (t, 4) = v n (t, sl A, for all t £ T as, 

Vn{t, 4 ) - 

max 
6n(r t )e 
B n {T t y n -,A) 

where a n (b s ) is given in Lemma 5.1, and 


f T - 1 


^2 d n (,b k ) - 6 n @l(s n - s„) >, (34) 


. k=t 


b„(7I,4;^) = 


{((/ 3nk,dnk),k £ Tt) S.t. p nk = ~ fn{d n k) + 25 n (T n - A) 

« and J2k=t d nk - min {-4, 0n(s„ - 4)}}, 

when t> 0 

^Bn(A), when t = 0, 

(35) 

with B n {A) defined in (16). The terminal value function is 
defined as, 

v n {T,sl)±-6 n Ql {sn-s T n ) 2 . (36) 

Note that the set of bids of the n-th EV over the interval Tt 
specified in (35) is defined in such a way that the total bidding 
demand over the whole interval T is guaranteed to equal A. 
The value function definition implies: 

t^n(0) s n , A, b_ n ) — max u n {b rLl b—n')i (37) 

bn£B„(A) 

and therefore, 

u n {b* n ,b- n ) = max v n (0 ,s°;A,b- n ). (38) 

Ae[o,r n ] 

Let b* (Tt) = (W* k ,d* nk ),k £ Tt) = b^iTt, 4; A, b- n ) 

denote the best bid of the n-th EV solving the optimization 
problem (34) over the interval Tt with respect to A and 6_„, 
and let n* (t, sl A, b- n ) denote the total bidding demand 
over the interval Tt of the n-th EV subject to the best bid 

b n {Tt, sl A, b_n), 

T-l 

K(t, si, A, 6_n) = Y, d nk(Tt, sl A, b-n). 
k=t 

Define B n {t, 4) = B n {t, sl A, 6_ n ), for any t £ T, as the 
set of bids at time t, such that 

B n {t, 4; A, b- n ) = 

|(/4i, dnt ) s.t. fj nt = —fn{dnt) + 2<5 n (T,, — A) , 
and d nt + n* (t + 1, 4 +1 ; A b~n) 

< min {A, 0„(s n — s^)}|, when t > 0 

|(/3 nt,d n t) S.t. /3 nt = — fn(d n t ) + 25 n (r„ — A) , 

and d nt + Hn(t + 1, sl 1 ; A,b^ n ) = a}, 

when t = 0, 

(39) 

where 4 +1 is given by (32). As with the set of bids of the 
n-th EV over the interval Tt specified in (35), the set of bids 
of at each time t £ T specified in (39) is defined such that the 
total bidding demand over the whole interval T is guaranteed 
to equal A. 

Theorem 5.1: The value function v n (t, sl A, b- n ) of the 
n-th EV, with respect to a fixed A £ [0,T n ] and a collection 
of bids b- n of the other EVs, can be implemented by solving 
the Bellman equation, 

Vn{t, S n ; A, b—n) = 

max i a n {b t ) + v n (t + 1, s^ +1 ; A, 6_ n ) \, 

bnt&Bn(.t,S^-,A,b- n ) l J 
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for t £ T, where B n (t, s^; A, b_ n ) is defined in (39), a n (b t ) 
is specified in Lemma 5.1, and v n (T,s x ) is given by (36). 

Proof. It is straightforward to verify Theorem 5.1 by apply- 
ing the optimality principle for the underlying optimization 
problem defined in (37) for the n-th EV with respect to A 
and 6„ n . ■ 

Finally, the n-th EV’s overall optimal bid, with respect to 
the bid profiles b n of the other EVs, is given by (38). 


B. Update mechanism for EVs 

The bid-profile update process of Section V-A motivates a 
heuristic algorithm whereby each EV successively updates its 
optimal bidding strategy with respect to the latest available 
bidding strategies of all the other EVs. If this iterative process 
converges, the resulting bid profile will be a NE, as established 
by Definition 1. 


Algorithm 1: (Nash equilibrium implementation algorithm) 

. Provide an initial bid profile b(°\ 

• Set the iterative step k = 0. 

. Set the iteration termination criterion e > e 0 for some 

e 0 > 0. 

« While e > e 0 

- For n = 1 : N 

Determine the best response for the n-th EV, 
b { n k+1) , with respect to b[ k+1) , b {k ^\ b { n % 

■■■ , bffl, by maximizing the payoff function, 

°n — argmax u n (o n , o 1 . , o n _ 1 ,o n+1 , ,o N ), 

bnGBn 


which can be achieved by applying the method 
developed in Section V-A. 

- Update e := | |fo( fc -t- 1 ) - fo( fc )|| 1 . 

- Update k := k + 1. 

End of Algorithm 

Whilst convergence of Algorithm 1 cannot currently be 
guaranteed, such guarantees have been established for related 
auction games. Various static iterative processes have been 
designed to determine the efficient NE for PSP auction games 
with a single type of resource, see for example [30], [31], 
[48]. Also, the so-called quantized-PSP auction mechanism 
was developed in [49]— [5 1 ] for electricity sharing games in 
a single time interval, with that process converging to a 
nearly efficient solution. Accordingly, ongoing research is 
studying the convergence properties of Algorithm 1 and the 
performance of the resulting NE for auction games where 
multiple resources are coupled. 


C. Numerical illustration 

To illustrate the auction-based coordination process, a nu- 
merical example will consider EV charging over a common 
time horizon T = 24 hours, from 12:00 on one day to 
12:00 the next day, with a time-step of AT = 1 h. The 
background demand for the example is shown in Fig. 4. For 
the purpose of demonstration, a small population of 5 vehicles 
will be considered. Each EV has a common battery capacity 
of 30 kWh and a common maximum SoC value = 0.9. 



Fig. 4. Background demand and the aggregate optimal charging strategies. 



20:00 0:00 

Charging Intervals 


- step 3 
step 4 

- step 5 

- converged strategy 

- basic demand 


Fig. 5. Convergent updates of Algorithm 1. 


Heterogeneity is introduced by letting the initial SoC values 
s (} n for the five EVs, prior to the charging interval, take the 
values s 0 = [0.1 0.15 0.23 0.14 0.08] T . 

The generation cost is given by Ct(x t ,D t ) = 
0.005 (Jfn.eN Xnt + -Dt) ‘ and the battery degradation 
cost by f n (x n t) = 0.002x‘^ lt . Both these functions are strictly 
convex. The weighting factor for the quadratic charging 
deviation cost of each EV is set to 5 n = 10 for all n € N . 

The efficient EV charging trajectory, given by Problem 1, 
is shown in aggregation (J2neN x nU t £ T) in Fig. 4. In 
contrast, the distributed approach to charging coordination, 
described by Algorithm 1 , gives the update evolution shown in 
Fig. 5. At each iteration, all EVs determine their optimal bid 
profile by solving the dynamic programming problem formal- 
ized by Theorem 5.1. Fig. 5 shows the aggregate allocation 
(^ ngA f t G T) obtained by the auctioneer solving 

Problem 2 with respect to the bid profile b'- kl at the fc-th 
iteration. Comparing Figs. 4 and 5, it is clear that the auction 
game converges to the efficient charging solution. 

At the efficient strategy b* , the aggregated payments of the 
EV population amount to XXeAf r? i(b*) = 23.44, and the 
increased generation cost is { Ct (-Dt+XXe N x lti b *))- 

c-t ( ff ) } = 22.48. The resulting surplus of 0.96, which is 4.1% 
of the total payments, can be used to reimburse the services 
provided by the system operator. 

VI. Conclusions and Ongoing Research 

As sales of electric vehicles (EVs) increase, the energy 
demanded by EV charging will begin to impact power system 
operation. Accommodating large numbers of EVs on the grid 
will require coordination of their charging requirements. A 
distributed approach to coordinating EV charging, based on 
the progressive second price (PSP) auction mechanism, has 
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been proposed. This auction mechanism ensures incentive 
compatibility. It is shown that the efficient (socially optimal) 
charging schedule is a Nash equilibrium (NE) of this auction. 

Individual EVs are capable of shifting their charging re- 
quirements in time across the charging horizon. Such cross 
elasticity complicates verification that the efficient solution is 
a NE. This has been addressed by partitioning players’ bids 
into subsets, each of which is composed of bids that possess 
the same total desired energy over the horizon. Consequently, 
cross elasticity is eliminated for bids within each subset. With 
this construction, the NE property of the efficient solution has 
been established by verifying for each subset that no player 
can benefit by unilaterally deviating from their efficient bid 
profile. 

An iterative update mechanism, based on dynamic program- 
ming, has been developed for the underlying auction game. 
This update mechanism was illustrated using a numerical 
example which showed convergence of the auction to an 
efficient NE. 

The developments in this paper motivate various ongoing 
research directions. This current work shows that for a single- 
step time horizon, the efficient solution is the unique NE for 
the underlying game. It remains to verify whether efficiency 
loss may occur for multi-step time horizons due to non- 
uniqueness of the NE. Work is also required to establish 
convergence properties of the auction update algorithm. 

The model underpinning the charging coordination scheme 
assumes that the EV population and the background demand 
are known with certainty prior to the charging period. The 
resulting auction game is static and the auction process can 
be undertaken off-line ahead of actual charging. In reality, the 
required information may be difficult to predict accurately, 
for example EVs may come and go without any advanced 
warning. Consequently, the efficient NE determined ahead of 
the charging interval may be suboptimal relative to a solution 
that considers disturbances. A practical approach to addressing 
this deficiency is to adopt a receding horizon strategy [52]. 
This is the focus of ongoing research. 

Appendix 

A. Proof of Lemma 3.2 

The bid price /3* t = g^—^w n (x**) can be written, 

fiu = 2 S n (T n -' 52 x%)-&(x%). 
teT 

Because P* t > 0 and f' n (x*ff) > 0 by Assumption A2, then 
J2teT x nt -Tn < 0. Together with (J2 tGT x nt -T n ) A n = 0 
in (5b), this gives A** = 0. 

With /3* t = gj—w n (x^*) in (5c), the complementarity 
condition (5a) gives, 

= C t ( D t + X nt ) . when X nt > 0 (40) 

< C t (D t + ff n eN x n*t ) ) when x * n * t = 0 , 

for all n £ Jf 1 1 e T. Substituting d* t = xf* t into (8b) gives, 
a* t > 0, « r <>; t = 0. (41) 


If x * \ = 0 then x* n t = 0 = x * : * t . If x** t > 0 then (40) indicates 
that /3* t = c' t (D t + Y,neAT x nt)- T he following argument 
shows that x* t < x * * t cannot occur. 

By the convexity property of Ct(-) and if x* t < x ** t , then, 

c’ t (D t + £ x* t ) - c' t (D t + <t) < 0- (42) 

nEA/* nEA/" 

Also, if x* t < x*j then (41) indicates that a* t = 0. Together 
with > 0 in (8), this gives, 

c[(D t + ]T x* nt ) - c' t (D t + x *nt) > 0, 

nEA/" nE M 

which contradicts (42). Hence, x* t < x** t cannot occur, 
implying x* t = x*) . In summary, under the bid profile b*, the 
associated optimal allocation x*(b*) for Problem 2 is identical 
to the socially optimal solution to Problem 1. 

Also, (10) follows from (40) with x* = x**. End of Proof 

B. Proof of Lemma 4.3 

If d* t = 0 then d n t = 0, according to (22). In that case, 
x n t = x* t = 0 trivially. The remainder of the proof assumes 

d* nt >0- 

The first step is to establish u n (b n ,bf n ). From (18), 

dnt = Put (dnt ; A) > Pnt(dnf,^2 d nt)> 

ter 

because A < Al so > the constraint d n t < d] j t 

ensures that, 

Pnt(dnt ; d nt ) > Pnt(d nt \ d nt ) = fi nt . 

teT teT 

Therefore f r ,t > Put* as illustrated in Fig. 6. Also, from 
Lemma 3.2, /3* t > /3^ t for all m £ N \ {n}, so p nt > P* nt . 



Fig. 6. Relationship between /5 n t and when 0 < A < 

Because X n t ^ dnt d n ^ — %nt Xrnt — X mt 

for all m G J\f \ {n}, then, 

c' t (D t + J2 *w) < c t(D t + E *«) = (43) 

keAf keAf 

for all j £ Af where d* t > 0. From (8a),(8c), c' t (D t + 
HkeAT X kt)-Pmt+Pmt > 0. Therefore (43) implies a m t > 0, 
and hence from (8b) that x m t = d* nt . Likewise, because 
Pnt > Pnt . c.'t(Dt + YlkeAT x kt) - Pnt + Pnt > 0 implies 
a n t > 0, and therefore that x n t = d n t ■ Consequently, all 
players are fully allocated, and (24) is established. 
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The payment of the n-th player is given by (13), with the 
above argument indicating that x^ff n = x^ t for all m £ 
AC \ {n}. Therefore, 

T n (b n ,b*_ n ) = 

E { c t( D t + E dmt + d nt )-c t (D t + E dmt) }• 

£G T m^n m^n 

The payoff of player n then follows from (14) as, 

tln(bm n) “ Xn(b n ^ b_ n ). 

The auctioneer needs to find an optimal bid b* that maximizes 
the payoff u n (b n: b*_ n ). This is achieved by the following 
restatement of the optimization problem (23): 

Problem 3: 

max u n (b n ,b*L n ), 

b ri 

such that (22) is satisfied together with, 

'y ( d n t = A. (44) 

ter 

The Lagrangian for Problem 3 is given by, 

L s {d ni X,n) = u n (b n ,b*_ n ) + ^2 X t {dnt ~ d* n 2j 

ter 

T P ^ y d n t ) 
ter 

where X t and p are the Lagrangian multipliers associated with 
the constraints (22) and (44), respectively. 

The KKT conditions for Problem 3 are given by, 

d d 

— ~ L (d n , A, /i) ^ 0, d n t ^ 0, — L (d ni A, fi)d n t = 0, 

ddnt ddnt 

dnt d n t E ^t ^ 0, ( d n t — 0? 

^2 dnt — A = 0, 
ter 

for all t G T, and, 

—~—L s {d n , A, p ) = — fn(d nt ) — c’ t (D t + ^ d* mt + d nt ) 
ddnt m ^ n 

+ Xt + ft. (45) 

Since d nt — d^ t < 0, then = 0, and so, 

dntKt) = fnKt ) + c[(D t +Y / d*mt+ <() > p, 

m^n 

where g n t(d n t) is defined in (25). If d* t > 0 then (45) equals 
zero, and so g' nt (d* lt ) = p. Therefore (26) holds. 

Under Assumptions (Al, A2), the necessary KKT condi- 
tions are also sufficient for optimality. Consequently, (26) is 
sufficient for the optimal bid 6* of player n. 

End of proof. 


C. Verification of inequality (28) in Theorem 4.2 

For notational simplicity, x and x denote the optimal allo- 
cations with respect to (b n ,b*_ n ) and (b a . b*_ n ), respectively. 

Consider the allocations of all EVs at time t\ with respect 
to the bid profile ( b n ,b*_ n ) tl . Because d ntl < d* tl < d* tl < 
d? ntl , (18a) gives, 

Pnt-i = fint\ (dnt! i A ) > finti ( d ntl A ) = /3 ntl . 

Also, by Lemma 3.2, /3* t > /3^ t for all m £ AT \ {n} when 
d* t > 0. Therefore, /3 ntl > Using an argument similar 
to that following (43), it is straightforward to show, 

Xnt! = dnt 1 , X mtl = d* mtl for all m £ AT \ {n}. (46) 

Hence, at time t \ , all EVs are fully allocated with respect to 
( b n ,b *_ n ) tl . Similarly, with respect to ( b n ,b *_ n ) tl , 

Xnti = dnt 1 , X mtl = d* mtl for all to G AC \ {n}. (47) 

By (46) and (47), the difference in the payments of the n-th 
EV at time t\ with respect to (b ni b*_ n ) tl and ( b n ,bf n ) tl is 
given by, 

Ar„ tl A Tnti (( b n , b*_ n ) t 2j - r ntl ((b n , b*_J t 2j 
= c t (D t i + ^2 d*mt i + dnti ) 

m^n 

— Ct(D tl + ^2 d*mti + dnt i)- (48) 

m^n 

For the n-th EV at time !■>. the difference in payments with 
respect to ( b n ,b*_ n ) t2 and (6 n ,6!_ n ) t2 is given by, 

= 7 nt2 b_ n )t 2 ^ /7 nt2 n)^2^ 

— Ct(Dt2 ^ ^ 3Cmt2 ^nt2 ) 

m^n 

— ct(D t , 2 + E %mt2 %nt2 ) 

m^n 

+ y ^ Pmt2 ( x mt 2 ~ £771(2)- (49) 

m^n 

The last term of (49) can be simplified by recalling from 
Lemma 3.2 that all EVs, k £ AC, with d* kt2 > 0 share the same 
value for /3 kt2 - Denoting that common value by 6* t2 allows 
the last term to be expressed as ^t 2 'Em? n (Xmt 2 £777(2)- 

It follows from (27d) that for the 71-th EV, they difference in 
payments at times t f t\,t 2 , with respect to (b n ,b*_ n )t and 
(b n ,b*_ n ) t is, 

A T nt = T nt ((bn,b*_ n ) t ) - T nt ( (b n , b *_ n ) t ) =0, Vf 

(50) 

Thus, by (48)-(50), the difference in the payments of the n-th 
EV with respect to (b n ,bf n ) and (b n ,b*_ n ) satisfies, 

A T n = T n (b n , b*_ n ) - T n (b n ,b*_ n ) = A T ntl + A Tnt 2 . (51) 

The difference in utility of the n-th EV, with respect to 

(b n ,b*_ n ) and (b n ,b*_ n ), is given by, 

Aw n A W n (x n ) - W n (x n ) 
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— ^ n ( ^ ' Xnt T 5 n ^ ^ ^ X n ± 

teT teT 

V fnidnt) fn(dnt) V fn(,Xnt) fn(,Xnt 2 )- 

(52) 

By (51) and (52), the difference in the payoff of the n-th EV, 
subject to {b n ,b*_ n ) and (b n ,b*L n ), becomes, 

A u n = u n (b n ,b*_ n ) - u n (b n , b *_ n ) = Aw n - A r n . (53) 

To establish (28), firstly the case with d n t 2 = d^ t2 will be 
addressed, then the three cases d n t 21 d n t 2 G 72.,;, * = 1,2,3 
will be considered separately. 

Case I, d n t 2 d n t 2 <C d n t 2 d n t 2 

Because d nt2 = d * lt2 , 

Pnt 2 = Pnt 2 (d nt2 , A) > /3 n t 2 (d nt2 , 'y^ d nt ) = (3 nt2 . 

teT 

Likewise, with d* nt2 < d nt2 < d nt2 = d* nt2 , 

)nt 2 = Pnt 2 ( d n t 2 , A) > )nt 2 (^nt 2 * = @nt 2 • 

A similar argument to that used to establish (46), (47) for t\ 
shows that all EVs are fully allocated at t 2 with respect to 
both (b n ,b*_ n ) t2 and ( b n ,b*_ n ) t2 : 

Xnt 2 d nt2 d n i 2 , r m £ 2 — d rn t 2 for all ui G J\f \ {A} , 

Xnt 2 d n t 2 , %mt 2 — dmt 2 for all ui G TV" \ {A}. 

Substituting these allocations into (53) gives, 

A U n — fn(dnt) fn(dnti') T fn(dnt 2 ) fn(dnt) 

- (ct(D tl + ^2 d* mtl + d ntl ) 

m^n 

- Ct(D tl + ^2 dmt! + dntl) 

m^n 

+ C t (D t2 + ^ 2 d*mt 2 + dnt 2 ) 

m^n 

- C t {D t2 + ^ 2 d*mt 2 + dnt 2 2j 

m^n 

= 9nt\[dnt\) dntii^nti) H” Qnt 2 {^nt 2 ) 9nt 2 ( d nt 2 ) 

2 dntiidnt^ydnti ~ d n t2 + g n t 2 (d n t 2 ){dnt 2 ~ d n t 2 ) 
— ft^dnti d n ti T d n t 2 d n t 2 ) 

= 0 , 

where the inequality holds due to the convexity of g n t{) and 
the subsequent equality follows from (26). Therefore, (28) is 
satisfied in this case. 

Case II, d n t 2 , d n t 2 G 72 | 

The initial step in showing (28) is to determine the allo- 
cations of all EVs at time t 2 with respect to the bid profile 
(b n ,b*_ n )t 2 . Firstly, consider d nt2 G lnt(72i). Then, 

Pm 2 > Pht 2 = c' t (D t2 + ^2 d* mt2 + d)t 2 ) 

m^n 


> c't( D t 2 + E d m t 2 + d n t 2 ) , 

m^n 

so it follows from the KKT conditions (8) that x n t 2 = d n t 2 . 

Now consider the case with d n t 2 = d„ t2 , the upper bound- 
ary of 72i. In this^case, p n t 2 = <4(A 2 + Em#n d *mt 2 +d n t 2 ). 
Assume x n t 2 < d n t 2 . Then due to the convexity of Ct(-), 

c t{ddt 2 ~\~ y ' d m f 2 -\-x n t 2 ) < c t (D t2 + 'y ' d mt2 -\-d n t 2 ) = )nt 2 - 

m^n m^n 

But (8) then implies a n t 2 > 0 and therefore that x n t 2 = d n t 2 . 
Hence a contradiction, so x n t 2 = d n t 2 ■ 

If d n t 2 > d kt 2 fh 611 it can l> e shown by contradiction 

that E mi tn x™t 2 = 0. Assuming *™,t 2 > 0 gives, 

c t (^*2 ^ ] x kt 2 + x n t 2 ) > c t (Dt 2 + y ^ d kt2 ) > (3 mt2 , 

k^n k£j\f 

for all m G J\f \ {n}. But (8) then implies x m t 2 = 0 
for all m G N \ {n}, hence a contradiction. Alternatively, 
if d nt2 < J2keMdlt 2 th en it can b e shown, once again 
by contradiction, that Y^keAf%kt 2 = Y^keNd*kt 2 - Consider 
E/cea/* *^ fci 2 ^ Eke*d k t 2 - Then c' t {D t2 + J)keX^kt 2 ) > 
/3^ t2 for to G Af \ {n}, with (8) implying x m t 2 = 0, 
hence a contradiction. If YlkeM ^ kt 2 < J2keM d kt 2 ’ then 
c' t (D t2 + EfceAA^) < /3mt 2 > with (8) implying x„ lt2 = 
^mt 2 - l ea( ^ s t0 another contradiction, as = 

d* 


d m t 2 + dnt 2 > J2keAf d kt 2 - Summarizing, 


%nt 2 — d , nt 2 i E %mt 2 0 ? E %mt 2 

m^n m^n 

if d n t 2 < y ^ d k t 2 , 

keX 

Xnt 2 — d n t 2 , y ^ x rn t 2 — 0 , y ^ 

m^n m^n 

if d nt2 > y ] d kt2 . 

keX 


■ d n t 2 — y ^ d kt2 , 

keX 

(54a) 

■ d n t 2 > y ' d kt2 , 

keX 

(54b) 


Similarly, the above analysis also holds for the bid profile 

{b n ,b*_ n ) t2 so, 


dnt. 2 — d kt2 , 
keM 

(55a) 

d nt 2 > y~~^ d kt2 , 

keM 

(55b) 


^nt 2 — d nt2 , E %mt 2 ^ 0? E ^mt 2 

m^n m^n 

if d n t 2 < E ^ 

keM" 

%nt 2 d n t 2 , E %mt 2 — O 5 E ^mt 2 

m^n m^n 

if d n t 2 > y ] ^fct 2 ■ 
keM 

Substituting into (51) gives, 

Ar n = Ct(Dt ! + ^ ^ dmtj + dnti) ~ Ct{Dti + ^ ' d mtl + dnt) 

m^n rriy^n 

+ Ct(Dt 2 ~\~ ^ ^ %mt 2 + dnt^ + fiot 2 ^ ^ (%mt 2 Xmt 2 ^) 
rriy^n m^n 

- c t (Dt 2 + ^ 3'mt 2 + dnt 2 ) • 
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Because x n ti -4 x n t 2 — x n ti 4- w nt 2 au d x n t — w n t for all 
t ^ ti,t 2 , it follows that Yt x nt = Yt^nt, and so (52) 
becomes, 

AWn = fn(d n t 1 ') “1“ fn(,d n t 2 ) fn(dnt 1 ') fn{d n t 2 ) • 

Three subcases must be considered, depending on the rela- 
tive values of d nt2 , d^ and Yk^N d lt 2 - 

Case II. 1, d n t 2 < d n t 2 < Yk=i d kt 2 : I n this case, A u n 
defined in (53) is established using (54a) and (55a), 

A U n — (Jnti ( dnti ) Qnt 1 (dnti ) A fn (d n t 2 ) fn ( d n t 2 ) 

- Pot 2 (dnt 2 - d nt2 ) (56a) 

> 9nti (dnti )(dnt 1 ~ d n t j) 

+ (fnKt 2 ) + (dnt 2 - d nt2 ) (56b) 

= 9nti(d ntl )(d n t 1 — d n t 2 ) + g n t 2 (d n t 2 )(d n t 2 — d n t 2 ), 

(56c) 

where (56a) holds by the specification of g n .t{‘) given in (25) 
and substitution from (54a) and (55a); (56b) holds by the 
convexity of g n t{-) together with (27a), and the convexity of 
/„(•) together with (29); and (56c) holds by (10) in Lemma 3.2 
and (25). 

From (22), d* nt2 < d* nt2 , so g' nt . 2 {d* nt2 ) > {/ nt2 {d* nt2 ) due to 
the convexity of g n t 2 { ')• By construction, d* ntl > 0, so (26) 
gives g'nfy{d* nt . 2 )^> g' ntl (d * ntl ) = g. Therefore, because (29) 
ensures d n t 2 > d nt2 , (56c) gives, 

A u n fj,(d n t 1 d n t 1 4- d n t 2 d n t 2 ) — 0, (57) 

where the final equality holds by (27c). 

Case 11.2, d nt2 < YkeM d kt 2 — d n t 2 : In this case, Au n is 
governed by (54b) and (55a), giving, 

A U n — 9nti(d n ti) 9nt\{dntr) T fn(d n t 2 ) fn(d n t 2 ) 

Ct(C)t 2 4~ ^ ^ Xmt 2 T d n t 2 ) c t (D t2 T d n t 2 ) 

m^n 

+ p; t2 ^ x rnt2 (58a) 

m^n 

^ gnt 1 (d n t 1 )(d n t 1 — d n t i) + f n(d n t 2 ) (,d n t 2 — d n t 2 ) 

+ c' t (D t2 +Y / dlJ(d nt2 -Y, %mt2 d n t 2 ) 

/cGA/" m^n 

+ Ptt 2 E *mt 2 (58b) 

m^n 

= gnt 1 (d n t 1 ){d n t 1 — d n tf) + f n (d nt2 )(d n t 2 — d n t 2 ) 

4~ Pot 2 {dnt 2 ~ d n t 2 ) , (58c) 

where (58a) holds by (53) and (25); (58b) holds by the 
convexity of g n t(-) together with (27a), the convexity of f n (-) 
together with (29), and the convexity of ct(-) using (55a); 
and (58c) holds by (10). Proceeding as in (56c), (57) yields 
A u n >0. ^ 

Case II. 3, y ', ucz \ r dj.,_ 3k d n t 2 5 d n t 2 - In this case, A u n 
uses (54b) and (55b) to give, 

A U n — 9nti ( d n ti ) 9nti ( d n ti ) 4“ fn ( d n t 2 ) fn (d n t 2 ) 

— Ct{Dt 2 4- d n t 2 ) 4- Ct(D t2 + d n t 2 ) (59a) 


> 9nti (.d ntl ) {dnt-t d n t 1 ) 4“ f n {d nt2 )(d n t 2 d n t 2 ) 

+ c t(A 2 4- dnt2)(dnt2 d n t2 ) (59b) 

— 9nti (^nti ) {dnti ~ dnti) 

n(.d n t 2 ) + c'(D t2 +Y / d* mt2 +dnt 2 )) 

m^n 

x ( d n t 2 d n t 2 ) (59c) 

= gnt 1 (d ntl ){d n t 1 — dnti) 4" 9nt 2 {d n t 2 )(dnt 2 — d n t 2 ) 1 

(59d) 

where (59a) holds by (53) and (25); (59b) holds by the 
convexity of together with (27a), and the convexity of 

/„(•) and ct(-) together with (29); (59c) holds by the convexity 
of c t (-) with d nt2 > YkeJ£ d*kt 2 ; and (59d) holds by (25)- 
Proceeding as in (57) yields A u n > 0. 

Hence, A u n > 0 whenever d n t 2 , d n t 2 £ . 

Case III, d nt2 ,d nt2 £ Int(TZ 2 ) 

The situation where d n t 2 £ IZ 2 will be considered as two 
separate cases. Case III, presented here, discusses d n t 2 £ 
lnt(IZ 2 ), while Case IV addresses d n t 2 = d^t 2 ’ tbe upper 
boundary of 1Z 2 . 

Consider the allocations of all EVs at time t 2 with respect 
to the bid profile (b ni b*_ n )t 2 . If d n t 2 < Yk^N d lt^ then 
the argument presented in Case II can again be used to 
show that YkexXkt 2 = YkeM d * kt2 - Also, because 0„t 2 > 
Pnt 2 = 4(A 2 +E fce M d * kt2 ), (8) implies x nt , 2 = cL 2 . Similar 
outcomes hold for the bid profile (6„ , b*_ n )t 2 a ® d n t 2 < d n t 2 . 
Therefore, (54a) and (55a) are again applicable. 

However, if d nt2 > YkeM d* k t 2 ’ then because 0 n t 2 < 
Pni 2 = c' t (Dt + E^„C 2 + 4t 2 )> there is no guarantee 
that x n t 2 = d n t 2 . Whether or not (54b) holds depends on the 
comparison between 0 n t 2 and c' t (D t2 + Ym^n %rnt 2 + d nt2 ). 
Similarly, for the bid profile {b n t 2 , bf n t2 ), there is no guar- 
antee that (55b) holds. 

Three subcases must be considered for x t2 and x t2 , de- 
pending on the relative values of d nt2 , d nt2 and YkeAf d k t 2 - 

Case III.l, d nt2 < d n t 2 < YkeM d *kt 2 : Analysis of A u n 
in this case is identical to that of Case II. 1, so A u n > 0. 

Case III. 2, d n t 2 ^f*k£Af d kt 2 — dnt 2 • Because 0nt 2 
P*nt 2 , satisfying (8) for the n-th EV results in Yk£Af%kt 2 > 
YkeN d kt 2 ’ w ith equality holding only if x n t 2 = d nt2 = 
EfcGA/* d kt 2 and Ymjin x m t 2 = 0. If the inequality is strict, 
then c' t (D t2 + EfceW ) > Pfnt 2 for all m £ M\ {?r}, with 
(8) implying Xmt 2 = 0. Hence, 

^ ^ d^t 2 — %nt2 — d n t 2 i ^ ^ %mt2 = 

k£J\f m^n 

The applicability of (54b) reverts to a comparison between 

0nt 2 and c' t {D t2 + dnt2 ) • 

• If 0 n t 2 > c’t (Dt 2 + d n t 2 ) then it can be verified that (54b) 
holds. Thus, A u n > 0, since the analysis in this case is 
identical to that developed in Case II. 2. 
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• If (3 nt2 < c' t (D t2 + d n t 2 ), (54b) does not hold. Rather, 
A u n can be established using (53), (25), (55a) and (60), 

A U n — 9nt\ (dnti ) 9nt\ (^nti ) 5 n ^ ^ ^ %nt 

ter 

4“ &n ^ ^ ^ Xnt 4“ fn(Xnt 2 ) fn(dnt 2 ) 

ter 

Ct (79^ 2 4“ ^ ^ •Emt2 4" d rh f 2 ) A Ct(Dt 2 A X n tfj 

m^n 

A fio t 2 ^ l 

m^n 

5nti (^nti ){dnti — d n tf) + f n (d n t 2 ){ x nt 2 — d n t 2 ) 

4” Srat 2 (^ni 2 ){dnti ~ ^nii 4“ <7ni 2 — X n t 2 ) 

+ $>t 2 E Xm ‘ 2 ■*" C t(-^2 + E ^t 2 4" d ni 2 ) 
m^n m^n 

X ( ^nt2 ^ ^ %mt2 dnt2 ) (613.) 

m^n 

= 9nt i (^nti ) (^nti — d n t 1 ) 4" <7 n t 2 (d n t 2 ) (<7nti — d n t 1 ) 

(61b) 

A p(^d n t 1 t^nti) 4“ /^(<7nti dnt -^ ) (61c) 

= 0, 

where (61a) holds by the convexity of g nt (-) together 
with (27a), the convexity of /„(•) together with (29) 
and (60), the convexity of c t (-) together with (55a) and 
(60), and the concavity of — ^(Eter ‘ Crlt ~ r n ) 2 together 
with Lemma A. 1 specified below, recalling that b n , b,, £ 
S n (A) with XXer d «* = < Et^r^t- and that 

X^££=7"' %nt X^ter ^ n t — d n t 1 ~\-d n t 2 {d n t 1 A2t n i 2 ) A 0, 
(61b) holds by (10) in Lemma 3.2 together with (25); and 
(61c) follows the same justification as (57) though using 
(27a). 

Lemma A.l: Consider an allocation x n (b) = (x nt , t £ T) 
with respect to a bid profile b, such that d n t < 

Eter<f Then: 

-Tn) ^ > 9nt(d* nt ) > P, (62) 

for all t £ T, where g nt is defined in Lemma 4.3. 

Proof of Lemma A.L 

d 


dx r , 


dn f ^ ^ x n t j — 2J n ^F n y ' x n t^ 

t&T ' teT 

>26 n (r n ~J2 d: 


ter 

c! { J* 


rrif^n 


= g'ntKt) 

> Pi 


Lemma 3.2; (63d) holds by the specification of g nt in (25); 
and (63e) holds by Lemma 4.3 and the convexity of g nt . 

^ End of proof of Lemma A.l. 
Case III. 3, Y^k^H d *kt 2 — ^«t 2 < r Et 2 •' Using the same 
argument as in Case ITT. 2 gives x m t 2 = E,n# n ^mt 2 = 

0. Analysis of Au ra depends on the relative values of d n t 2 , 
x n t 2 and d n t 2 , keeping in mind from Lemma 4.1 that 

(5nt 2 (dnt 2 i 2^) A Pnt 2 (dnt 2 i 4l) - 

• If x nt2 = d nt2 then x nt2 = d nt2 must also hold. Analysis 
of Au r , in this case is identical to that developed in 
Ca£e II. 3. 

• If d n t 2 A x n t 2 < d n t 2 then x n t 2 — d n t 2 A x n t 2 . Analysis 
of Au n follows that of Case III. 2. 

• If x n t 2 r d n t 2 r d n t 2 then x n t 2 ^ x r it 2 , and A u n 
satisfies, 

P^ u n > 9nti (d n t 2 ) (dnti — 4" f n {Xnt 2 ){ x nt 2 — X n t 2 ) 

4~ ^t. (U^ 2 A x n tf){x n t 2 x n tfj 

4 “ 9nt 2 (. d ’nt 2 )(dnt 1 ~ d nt\ ) 

4“ 25 n ^r n y ^ x n t ^ (x n t 2 x n t 2 ) (64) 

t£T 

A /n^ra^ ) X n t 2 ) 

4“ 4~ x n t 2 )(x n t 2 x n t 2 ') 

4~ 2<5 n ^L n y ' Xnt^J {x n t 2 x n t 2 ) , (65) 

ter 

where (64) holds by the convexity of g ntl (•) together with 
(27a), the convexity of f n (-) and Ct(-), the concavity of 
— <5„(r n — Eter^t) 2 and Lemma A.l; and (65) makes 
use of (61b). 

Further analysis uses x n t < d n t for all t £ T to give, 


25 n (r n y ^ x n t J A 

teT 


) >2 j n (: 


r?9. ^ ^ d n t ^ 

ter 


— fn(dnt 2 ) 4“ Pnt 2 i 


( 66 ) 


(63a) 


= /?:* + /;«t) (63b) 

= <4(A + E + O + /™Kt) ( 63c ) 


(63d) 

(63e) 


where (63a) holds because yf ter x nt < X^ter^"* < 

E ter d nt5 (63b ) follows from = ah w n(d* n ) = 
-fnKt) + 2bn(r n - Eter^nt); (63c) holds by (10) in 


where the equality follows from (11). Because x n t 2 > 0 

and Em#n^™t 2 = ( § ) S iveS Pnt 2 > c' t (D t2 + X n t 2 ). 

Therefore, 

2‘5 n ^T' n y ^ d n t) A /n(^7it 2 ) 4~ A ^ n ^ 2 ) 

ter 

> fn(Xnt 2 ) + c' t (D t2 + J nt2 ). (67) 

Because ir„ t2 < x nt2 , (65) and (67) ensure A u n > 0. 
Hence, A u n > 0 whenever d n t 2 ,d n t 2 £ Int(72.2)- 

Case IV, d nt2 — d n t 2 , d n ± 2 £ 77.2 

In this case, fi nt2 = /3* t2 , so (8) ensures that EfceAf^s = 

EfceA/ - d kt 2 and d *nt 2 < Xnt 2 < &nt 2 - 

Case IV. 1, d nt2 < EfceAf^fet 2 ; Using the same argument 
as in Case III, (55a) j.s again applicable. 

If x nt2 > x nt2 = d n t 2 , dAu n can be established by, 

Au ra > g n t 1 {d nti )(d n t 1 — d n t 1 ) + f n {d nt2 )(x n t 2 — d n t 2 ) 

4" 9nt 2 ( d n t 2 ) (dnti ~~ d n t i + d n t 2 — x n t 2 ) 

4~ Pot 2 (*® 71*2 d n t 2 ) 


(68a) 
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— 9nti i^nti ){dnt 1 d n t 1 ) + g nt2 (d nt2 ) {d n t 1 d n t 1 ) 

(68b) 

> 0, (68c) 

where (68a) holds by the convexity of g n t{-) together with 
(27a), the convexity of /„(•) together with x n t 2 > x n t 2 
and (55a), and the concavity of —5 n (ff2 lt< ^-x n t — T„) 2 to- 
gether with Lemma A.l, recalling that b n ,b n £ B n {A) 
with Eterdnt = A < Eterd*,A (68b) holds by (10) in 
Lemma 3.2 together with (25); and (68c) follows from (61b). 
If x n ± 2 x n t 2 — d n t 2 , At/^ is given by, 

A u n > g n t 1 {d n t 1 ){d n t 1 — d n t i) + f n (d n t 2 )(x n t 2 — d n t 2 ) 

+ 9nt 2 (d n t 2 ){d n t 1 — dnti ) “b fiot 2 i^nt 2 d n t 2 ) 

“b 2b n (r n ^ ' %nt')(dnt 2 •^nt 2 ) (69a) 

tGT 

^ fn(dnt 2 )(,%nt 2 d n t 2 ) + ftot 2 ( x nt 2 d n t 2 ) 

~b 2b n (r n y ' %nt)(dnt 2 r n ± 2 ) (69b) 

ter 

> 0, (69c) 

where (69a) holds by the convexity of g nt ( •) together with 
(27a). the convexity of /„(■) together with (55a), and the 
concavity of — 4(EteT x nt~ r„) 2 together with Lemma A.l, 
recalling that b n ,b n £ B n (A) with Jf teT dnt = A < 
J2teT d nt’ (69b) uses (68b); and (69c) uses (66) together with 

Pnt 2 > /?o t 2 - __ 

Case IV.2, d nt2 > EkeU dtt 2 : Using tbe same argument 
as in Case III, (55b) is applicable. Then similar to the analysis 
of Case IV. 1, > 0. 


Case d n t 2 , d n t 2 £ 7^-3 

In this case, (3 nt2 < fd nt2 < /3* t2 ’ so (8) ensures that x m t 2 = 
x m t 2 = d* mt2 for all m £ Af\ {n}, and x nt2 < x nt2 < d* t2 . 4 
Hence, (51) becomes, 

A r„ = 

Ct{Dt i + Y, 4 t i + dnti ) - ct{D t i + Y, d* mtx + d ntl ) 

m^n m^n 

+ c t (D t2 + Y dmt 2 ~b x n t 2 ) ct(Dt 2 + 'y ' d mt 2 +r n t 2 ). 

m^n m^n 

(70) 


Using (52) and (70) in (53) gives, 

A ll n — Qnti(.d n t\) 9nti(dnti') “b gnt 2 (.Xnt 2 ') 9nt 2 {pCnt 2 ) 



ter 


+ d n 


{^2 Xnt 

ter 



> -ft 
- 6 . 


^d n ti d n t2j 9nt 2 (d n t 2 ) (%nt 2 %nt 2 ) 

i ^ ^ ' %nt “b S n ^ y ' X n t r r 


(71) 


4 The equality x n t 2 = x n t 2 = 0 can occur if c’ t (Dt 2 + J2 m ^n ( Knt 2 ) > 

0nt 2 > fnt 2 ■ 


where the inequality holds by the convexity of g nt (- ) together 
with (27a) for the first term, and with x nt2 < x n t 2 < d) \ t2 for 
the second term. 

Using (62) from Lemma A.l together with (27a), the 
concavity of S n {Y, teT x n t - r„) 2 , and x nt2 > x nt2 gives, 

-4 


^ %nt r n 

^ ~l~ fin f ^ v %nt 

ter 

K ter 

^ 9nt 2 (^nt2) 

^d n t 1 d n t 1 H- x n t 2 x n t 2 

[I (d n t 1 

dntl J 9nt 2 (d n t 2 ) ( x nt 2 ~ 


Thus, Jt follows from (71) and (72) that A u n > 0 whenever 

dnt 2 -)d n t 2 £ 7^.3. 

In summary, the analysis presented in Cases I-V shows that 
inequality (28) holds for all d nt2 > d* t2 . 

End of proof. 


D. Proof of Lemma 4.4 

The total aggregate payment made by all EVs is given by 

(13), 

Y r «( b *) = Y J2 Tnt ( b *')’ 

neN neUter 

with 

Tnt(b*) =-ct(D t +Y^ x ™t n ) + Ct(D t + X *kt ) 

m^n k£j\f 

+ Y^mt{ X ^ n - X *mt)}, (73) 

m^n 

where x* kf denotes the optimal charging allocation of the fc-th 
EV given by Problem 2 with respect to b ^ , and x* n f " denotes 
the optimal charging allocation of EV m ^ n with respect 
to (0 ni b*_ n ) t . It follows from the KKT conditions (8) that 
x mt n = X mt f° r ah rnf^n. Therefore, (73) simplifies to, 

T nt {b* t ) = C t (D t + Y X kt ) - C t(A + Y X *mt)- 
k£Af m^n 

Because 0 < '£ m ^n x *mt < T,keAf x kv Assumption (Al) 
implies, 


Ct{D t + EfceV X kt) ^ C (A + Em^u x *mt ) 


> 


Ct(A + Ek£Af X *kt) -c{Dt) 


E 


feeAf x kt 


Therefore, 

Y r ™( b *) 

n£j\f 


Ct {d2t + EfegAf x kt) Ct(Dt + x mt) ^ 

x Xnt 

■neuter nt 


^ ct(D t + EfcgA r x kt) ct(Dt) _ * 

> Z^ Z^ v V* Xnt 

new ter e^keU kt 


c t{Dt+ YjkeU x kt) Ct(D t ) ^ ^ 

- 2^ v T* E, Xnt 

2^keU x : 


ter 


keU * kt 


n.eM 


^2 \ct(D t + ^2 x *t) - ct(A) J. 


ter 


keU 


End of Proof. 
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E. Proof of Theorem 4.4 

It will be shown initially that for each EV, n £ Af, there 
always exists a best response b n = (/3 n ,d n ) under which it is 
fully allocated, x n = d n . Denote by &*(&_„) = (/3*,d*) the 
best response of the ro-th EV with respect to b - n , 

K( b -n) = argmax u n (b n , 

S n eB n 


and let x* be the allocation with respect to (b* . b _„)• As 
stated in Section III-A, incentive compatibility holds, so a bid 
with price satisfying /3 n = w' n (d n ) is the best choice among 
all possible bids. Thus, it will be assumed all bids are truthful. 

Consider another bid b n = (p n . d n ) with f3 n = 4(4) 
and d n — X n , and denote by x the allocation with respect to 
( b n ,b - n ). The concavity of w n (-) and x* < d* ensure that 

4(4) > w' n (d* n f which implies that P n >P*. 

Assume x n < d n = x„ . Then, 

c'{D + 'ffx* m + x n )< c\D + ^fx* m + x* n ) <P*< p n - 

m^n m^n 


But (8c) then implies a n > 0 and so x n = d n . Hence a 
contradiction so x n > d n . But x n < d n , so x n = d n = 
That is to say, the bid b n results in full allocation. Since b n 
has the same allocation as 6*, it follows that b n has the same 
payoff as 6*. Therefore b n is also a best response. It may be 
concluded that for each EV, there exists a best response with 
a full allocation. Hence, the remainder of the proof assumes 
that EVs consider their best response with full allocation. 

Suppose there exists another NE, denoted by fo°, which 
differs from b * . By Lemma 3.2, b° is inefficient and does 
not satisfy (10). Also, according to Definition 1, 6° is the best 
response of the n-th EV with respect to 5° n . From the above 
argument regarding full allocation, it may be assumed that 
x^ = d 1 ^ for all n £ JV. Moreover, since b° satisfies (8): 

oO f= C '( D + E fce v4) + a n, if 4 > 0 

Pn \<c'(D + j: keM x°) + a n , if 4 = 0, 

where o n > 0. Since 4 = 4 and 6° * s inefficient, at least 
one EV, l £ A f, must satisfiy ct; > 0. 

Consider another bid for the Z-th EV, denoted by bi = 
(, fii,di ), such that c'(D + X^fceJV 4) < Pi < Pi- L et x be the 
allocation with respect to {bi,b°_fj. From (11), the concavity 
of w n (-) gives d° < di . Then, because 4 j = f° r a H 
and d(D + J2kesS x k) < i 1 follows that: 

0 < 4 < 4 (75a) 


Because xi > 0, the KKT conditions (8) require: 

c'(D + y x k ) < Pi- 
keJT 


Furthermore, (8) implies: 


= x 

< X 


0 

m 5 
0 


if P° m > Pi, 

otherwise, 


for all to f l, 


and 


X] 4 > x °k- 

k£Af keAf 


(75b) 


(75c) 

(75d) 


Hence the difference of payoffs between b° and (bi, b°_ ; ), 
given by A ui = ui(b°) — ufb^b^ff), satisfies the following 
analysis: 


Au, = Mxi) - - <5* ((4 - r,) 2 - {xi - Ti) 2 ) 

+ c{ d +Y1 4) - c(D + 4) 

fceAf feGAT 

+ ^2^m{ X °m- X m) 
rrif^l 

< f[(xi)(xi - 4) + 2 Si{xi - T t )(xi - 4) 

+ c'(d + 4 ) ( 4 - 4 ) 


fceAf fceV fceV 


■?<E' 




m^l 


(76a) 


< (// (dj) + 25i{di - TO) (xi - 4) 

+ c'(d + 4) ( 4 - 4) 

feeV fceV feeV 

+ Pi ^ ] i. x m ~~ x m) 

m^l 

= (c'(d + 4) - a) ( 4 - 4) ( 76c ) 

keAf keAf keAf 

< o, (76d) 


(76b) 


where (76a) holds because of Assumptions (Al, A2) and 
(75c); (76b) follows from Assumption (A2), xi < di and (75a); 
(76c) uses (11); and finally (76d) uses (75b) and (75d). 

This analysis implies that upb 0 ) < upbi, bff) which 
contradicts the statement that Ip is the best response of the 
l-th EV with respect to bfj. Hence, no NE exist apart from 
b*. The NE for single-interval auction games is unique and 
efficient. End of proof 
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